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Abstract 


A  program  of  investigation  of  theoretical  equations  of  state,  and  ultra¬ 
sonic  measurements  to  10  kbar  and  800°K  is  described.  The  theoretical  work 
summarized  includes:  reduction  and  analysis  of  shock  wave  data  for  rocks  and 
minerals,  formulation  of  complete  p~V-T  equations  of  state  for  calculating 
high  pressure  and  high  temperature  Hugoniots,  release  adiabats,  and  isotherms, 
and  the  examination  of  the  empirical  relations  between  bulk  modulus  and 
compressional  velocity  and  density.  The  application  of  Kroner's  theory  of 
random  aggregates  to  the  prediction  of  rock  elastic  properties  using  single¬ 
crystal  elastic  constant  data  is  discussed.  Ultrasonic  elasticity  measurements 
on  single  crystal  spinel  (MgAl^O^)  at  ambient  conditions,  and  on  rock  salt  (NaCl) 
to  pressures  of  10  kbar  and  800°K  are  also  summarized. 

1.  Introduction 

Equations  of  state,  or  constitutive  equations,  provide  a  complete  mechanical 
as  well  as  thermodynamic  description  of  material.  An  equation  of  state  relates 
stress  or  pressure,  strain  or  density,  and  one  or  more  thermodynamic  variables, 
such  as  internal  energy,  or  temperature,  for  a  specific  material.  Equations 
of  state  for  earth  materials  are  necessary  for  the  calculation  of  intense 
stress-wave  propagation  and  seismic  coupling  resulting  from  explosions  or 
impacts  on,  or  within,  the  earth's  crust.  Because  of  the  wide  range  of  dynamic 
pressures  which  are  induced  within  the  vicinity  of  a  nuclear  explosion, 
the  equations  of  state  which  are  of  interest  will  describe  the  response  of 
a  medium  over  a  tange  of  conditions  from  thousands  of  kilobars  and  tens  of 
thousands  of  degrees  down  to  a  fraction  of  a  bar  and  ambient  temperature;  the 


1. 


latter  conditions  correspond  to  large  distances  from  the  disturbance. 

The  present  program  which  is  a  continuation  of  a  previous  program 
summarized  by  Ahrens  (1970)  has  been  theoretically  concerned  with  enuations 
of  state  of  earth  materials  under  the  wide  range  of  dynamic  pressures  and, 
to  a  lesser  degree,  with  gathering  some  ultrasonic  data  at  high  temperatures 
(800  K)  and  moderate  pressures  (10  kbar)  for  some  earth  materials. 

The  theoretical  phase  of  our  program  has  concentrated  in  the  following 
areas : 

(a)  Analysis  of  the  increasing  quantity  of  shock  wave  data  for  silicates 
and  silicate-bearing  rocks. 

(b)  Formulation  and  application  of  a  complete  and  thermodynamically 
consistent  p-V-T  equation  of  state  suitable  for  calculating  Hugoniots, 
adiabats.and  isotherms  at  very  high  compressions  and  temperatures.  This 
class  of  equations  of  state  are  based  on  finite  strain  theory  as  carried 
out  to  the  fourth-order  term  in  Helmholtz  free  energy. 

(c)  Formulation  of  systematic  relationships  between  the  bulk  modulus 
and  compressional  velocity  with  density  for  a  variety  of  minerals 
and  oxides. 

(d)  Preparation  of  several  topical  reviews  on:  equations  of  state  and 
strength  of  material  data  pertinent  to  seismic  coupling  and,  ultrasonic 
and  dynamic  and  static  pressure-volume  data  for  rocks  and  minerals. 

In  the  experimental  phase  of  our  program,  we  carried  out  high  pressure 
and  high  temperature  measurements  on  the  elastic  properties  of  rock  salt.  In 
addition, some  high  quality  ultrasonic  data  on  the  elastic  properties  of  the 
mineral  ,  spinel  were  obtained. 
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2.  Research  Results 

The' papers  which  resulted  from  research  fully  and  partially  supported 
under  this  contract  were : 

1.  Application  of  the  Fourth-Order  Anharmonic  Theory  to  Prediction  of 
equations  of  State  at  High  Compressions  and  Temperatures,  submitted  to 
Physics  of  the  Earth  and  Planetary  Interiors.  T.  J.  Ahrens  and  L.  Thomsen. 

2.  On  the  Effects  of  Pressure  upon  Rock  Elasticity,  submitted  to  Physics 
of  the  Earth  and  Planetary  Interiors.  L.  Thomsen. 

3.  Elastic  Properties  of  Minerals,  R.  C.  Liebermann  and  E.  Schreiber, 
Transactions,  American  Geophysical  Union.  52.  142-147,  1971, 

4.  Static  and  Dynamic  Compression  of  Earth  Materials,  T.  J.  Ahrens  and 

T.  Takahahsi,  Transactions,  American  Geophysical  Union.  52 .  147-153,  1971, 

5.  Recent  Progress  in  the  Study  of  Dynamic  Rock  Properties  Pertinent  to 
Predicting  Seismic  Coupling,  T.  J.  Ahrens,  submitted  to  Proceedings . 

1970  ARPA  Seismic  Coupling  Symposium . 

6.  On  the  Elasticity  of  Polycrystals  ari  Rocks,  L.  Thomsen,  J.  Geophvs.  Res., 
in  press,  1972. 

7.  The  Bulk  Modulus-Volume  Relationship  for  Oxides,  D.  L.  Anderson  and 
0.  L.  Anderson,  J.  Ceophys ■  Res..  75.  3494-3500,  1970. 

8.  Velocity-Density  Relations,  D.  L.  Anderson,  J.  Geophvs.  Res..  75, 
1623-1624,  1970. 

9.  Revised  Shock-Wave  Equations  of  State  of  High-Pressure  Phase  of  Rocks 
and  Minerals,  G.  F.  Davies  and  D.  L.  Anderson,  J.  Geophvs.  Res..  76, 
2617-2627,  1971. 

10.  Shock  Wave  Equations  of  State  of  Minerals,  T.  J.  Ahrens,  Proceedines. 
International  School  of  Physics.  Enrico  Fermi,  in  press,  1972. 
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11.  Equation  of  State  of  NaCl:  Ultrasonic  Measurements  to  8  kbar  and 
800°K  and  Static  Lattice  Theory,  H.  Spetzler,  C.  G.  Saimnis,  and  R.  J. 
O'Connell,  submitted  to  The  Journal  of  the  Physics  and  Chemistry  of 
Solids.  1972. 

12.  Quasi-Harmonic  Finite  Strain  Equation  of  State,  G.  F.  Davies,  to  be 
submitted  for  publication,  1972. 

The  significant  results  we  obtained  (and  reported  in  these  papers)  are 
summarized  in  the  following  discussions: 

A.  Formulation  of  Equations  of  State  and  Theoretical  Analysis  of 
Experimental  Data 

Using  a  revised  value  for  the  relation  between  the  zero-pressure  seismic 
parameter,  4>q  =  j  ,  density  p^,  and  mean  atomic  weight,  M,  given  by 

P„  (g/cm3)  =  0.0492  M  $1/3  (1) 

o  o 

where  p  is  pressure,  p  is  density  and  s  is  entropy,  Davies  and  Anderson 
(paper  2)  reduced  the  high-pressure  shock  wave  data  for  a  series  of  19  rocks 
and  minerals  to  adiabatic  conditions.  In  a  revision  of  an  earlier  paper 
(Ahrens  et_  al . ,  1969)  they  found  that  the  adiabats  could  br  fit  using  a 
Birch-Murnaghan  equation.  This  equation  predicts  a  smooth  increase  of  4> 
with  increasing  density  (or  pressure). 

The  resulting  adiabats  are  derived  from  the  shock  wave  data  using  a 

Mie-Gruneisen  equation  of  state.  In  equation  (1)  $  is  given  In  units  of 
2 

(km/sec)  .  For  three  iron-rich  minerals  the  emprical  equation 

PQ(g/cm3)  =  0.051  M$1/3  (2) 
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for  the  relation  of  density  and  seismic  parameter  at  zero  pressure  was 
applied.  The  resulting  p-p  equations  of  state  are  observed  to  be  roughly 
parallel  in  the  log  (p/M)  versus  log  4>  plane  (Figure  1). 

Upon  examining  the  relation  of  bulk  modulus  and  density  (actually 
specific  volume)  Anderson  and  Anderson  (paper  7)  found  that  the  product  of 
bulk  modulus  and  specific  volume  for  different  classes  of  oxides  is  found  to 
be  effectively  a  constant  for  a  ^iven  set  of  structures  and  coordination  numbers. 

A  second  useful  result  in  an  empirical  relation  has  been  found  between 
the  reduced  Madelung  constant,  and  the  ionic  packing  in  a  lattice.  This  result, 
when  applied  with  knowledge  of  the  bulk  modulus,  can  be  used  for  calculating 
theoretical  equations  of  state  in  materials  having  potentials  of  the  essentially 
ionic  type. 

A  large  volume  of  both  ultrasonic  and  shock  wave  data  describing  the 
elastic  constants  and  p-V-T  equation  of  state  have  been  reported  in  the  literature 
for  single  crystals  and  monomineralic  samples.  Thomsen  in  a  series  of  two 
papers  (2  and  6)  has  examined  the  application  of  Kroner's  (1967)  theory  to 
prediction  of  the  elastic  constants  and  equation  of  state  of  polycrystalline 
materials  and  rocks  (composed  of  more  than  one  mineral)  in  terms  of  the  elastic 
moduli  of  the  individual  crystal  species.  In  the  case  of  polymineralic  rocks, 
the  lack  of  complete  data  make  tests  of  Kroner's  perfectly  disordered  aggregate 
theory  versus  the  more  common,  but  empirical,  Reuss-Voigt-Hill  average,  difficult. 
However,  the  kroner  theory  provides  a  self-consistent  framework  for  extrapolation 
of  near-zero  pressure  data  to  high  pressure  (paper  2).  In  the  case  of  a  highly 
porous  aggregate  this  theory  appears  to  describe  the  available  data  well 
(Figure  2).  It  is  in  this  context  that  the  theory  will  probably  be  most  useful 
in  predicting  elastic  properties  of  porous  rocks  on  the  basis  of  their  composition 
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Fig.  1.  Log  4  versus  log  (p/M)  for  adiabats  of  17  high-pressure  shock- 
induced  phases.  (a)  zero-pressure  constraint  on  adiabats  given  by 
p  =  0.048  M  *0.323  as  previously  used  by  Ahrens  et  al .  (1969). 

Short  dashed  curves,  oxides;  long-dashed  curves,  iron-rich  compounds, 
dash-dot  curves,  rocks  containing  feldspars.  (b)  zero-pressure  con¬ 
straint  on  adiabats  given  by  Eqn.  2.  Symbols  for  rocks  and  minerals 
are  identified  in  Table  1  (after  Davies  and  Anderson,  1971,  paper  9). 
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Figure  2 


Normalized  modulus  A  +  2u  =  v  j.  c 

from  Warren  fiq  q  c  4u/3  for  porous  obsidian;  data 

rrom  Warren  (1969),  theory  after  Kroner  and  ^ashin-qitriknvm 
Agreement  with  Kroner’s  theorv  ic  f  ,  trikman. 

f  (after  Thomsen,  1972)  Y  £°r  shear  “dulus 
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and  density  alone. 

Recent  progress  in  the  study  of  dynamic  yielding  of  solid  and  porous  rocks, 
the  effect  of  water  on  the  equation  of  state  and  very  high  pressure  equations 
of  state  was  reviewed  by  Ahrens  (paper  5  )  at  the  1970  ARPA  Symposium.  New 
Hugoniot  data,  both  above  and  below  the  Hugoniot  elastic  limit  was  compared 
with  hydrostatic  compression  data.  These  comparisons  indicate  that 
appreciable  stress  differences,  comparable  to  those  existing  at  the  dynamic 
yield  point,  are  supportable  by  rocks,  such  as  sandstone,  limestone  and 
granite  above  the  dynamic  yield  point.  The  quasi-static  failure  tests  which 
have  recently  been  reported  provide  data  which  closely  satisfy  a  Prager-Drucker- 
type  yield  surface.  In  particular,  the  one-dimensional  stress  tests  for  a 
series  of  porous  and  non-porous  rocks  indicate  sensitivity  of  fracture  stress 
and  dynamic  (Young's)  modulus  on  strain  rate.  Of  special  importance  is  the 
observation  that  water  affects  the  dynamic  flow,  and  resulting  stress  waves 
from  underground  explosions  because  of  steam  formation.  Water  strongly  affects 
the  phenomenon  of  block  sliding  which  takes  place  after  the  initial  shock  has 
propagated  through  an  initially  competent  rock.  This  mechanism  is  believed 
to  represent  the  dominant  process  limiting  shear  stresses  at  late  times  in 
the  flow  around  explosions  in  originally  strong  rocks.  He  also  notes  that 
at  high  shock  pressures  a  series  of  shock-induced  phase  transitions  involving 
changes  of  from  10  to  60%  in  density  takes  place  in  silicate  minerals.  These 
phase  changes  account  for  most  of  the  compression  which  occurs  in  the  first 
megabar  of  pressure  in  solid  rocks.  As  a  result  of  these  transitions  the 
release  adiabats  upon  initial  unloading  from  Hugoinot  states  are  considerably 
steeper  in  the  stress-density  plane  than  the  corresponding  Rayleigh  lines  and 
give  rise  to  appreciable  intrinsic  shock  attenuation. 
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In  1970,  a  fourth-order  finite  strain  equation  of  state  which  self- 
consistently  treated  the  thermal  properties  of  solids  was  first  described  by 
Thomsen  (1970) .  One  important  merit  of  this  theory  is  that  the  Gruneisen 
parameter  (y)  which  is  so  vital  in  adequately  describing  high  pressure 
Hugoniots  and  release  adiabats  for  solids  is  explicitly  given  as  a  function  of 
volume  in  terms  of  other  thermodynamic  properties  and  the  value  of  y  at  ambient 
conditions.  The  theory,  originally  derived  in  order  to  describe  p-V-T  equations 
of  state^has  later  been  extended  so  as  to  describe  shear  moduli  by  Thomsen 
(1972).  A  major  disadvantage  of  the  finite  strain  approach  is  that  the 
definition  of  strain  is  not  unambiguous.  Originally  Thomsen  derived  an 
equation  of  state  by  expanding  the  Helmholtz  free-energy  in  terms  of  the 
Lagrangian  strain  parameter  n, 


n 


i  r  v  \  2/3  1 

ur;  -1! 


(3) 


where  V  is  the  volume  of  the  solid  in  the  rest  state  (at  0°K  with  the  effect 
of  zero-point  vibrational  energy  removed).  This  theory  was  briefly  summarized 
in  our  previous  Final  Report  (Ahrens,  1970).  More  recently  Davies  (paper  12) 
derived  another,  more  general,  and  probably  equivalent  form  of  the  equation  of 
state  using  strain  measures  such  as 
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where, rather  than  relative  to  the  rest  state,  the  standard  volume  is  taken, 
more  conveniently,  under  ambient  conditions.  Because  of  the  above  mentioned 
ambiguity  in  the  definition  of  the  finite  strain  the  various  finite  strain 
theories  begin  to  diverge  from  each  other  at  compressions  V/Vq  £  0.7.  In  spite 
of  this  drawback,  because  of  the  potential  usefulness  of  finite  strain  theory 
in  generating  theoretical  Hugoniots,  isotherms,  and  release  adiabats  for  a  wide 
class  of  both  initially  solid,  and  porous  materials  using  low-pressure  ultrasonic 
and  thermodynamic  data,  Ahrens  and  Thomsen  (paper  1)  exhaustively  tested  its 
applicability  for  17  metals,  compounds, and  minerals.  The  materials  studied 
were  all  cubic  (crystal  structure)  and  ones  for  which  good  elastic  constant, 
thermodynamic  and  shock  wave  data  were  available.  Ultrasonic  values  of  the 
bulk  modulus,  K,  its  pressure  derivative  K’  =  dK/dP  and  in  some  cases,  the 
second  pressure  derivative  of  bulk  modulus  K"  =  dK'/dP  was  employed  in  the 
analysis. 

They  concluded  that  the  fourth  order  anharraonic  theory  is  adequate  to 
predict  P(V)  for  the  metals  to  compressions  of  (AV/V^)  x  »  0.3,  and  in  a  few 
cases  to  much  greater  compression.  For  ionic  compounds,  the  maximum  com- 
oression  is  generally  smaller,  and  more  erratic  within  each  croup  of  halides. 

The  uncertainties  associated  with  the  two  covalent  compounds,  MgO  and  BeO, 
make  generalizations  for  this  class  hazardous.  Experimental  resolution  of  the 
experimental  uncertainties  in  the  parameters  K/  and  K^Kq  will,  eventually 
allow  a  more  definitive  test  of  the  fourth  order  theory  for  this  class  of 
materials.  The  purely  thermal  contributions  to  the  pressure  seem  to  be  well 
described  for  the  metals  up  to  their  melting  points  and  for  MgO  to  temperatures 
almost  three  times  the  Debye  temperature.  However,  the  description  of  the 
free  energy  function  as  a  function  of  temperature  can  be  improved  by  using  a 
more  refined  model  for  the  specific  heat  (Kieffer  and  Kamb ,  1971).  The  elastic 
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moduli  are  described  well  to  compressions  of  perhaps  0.2,  although  the  data  is 
very  limited. 

On  this  basis,  they  concluded  that  the  anharmonic  theory  is  probably 

adequate  to  describe  the  elasticity  of  all  cubic  minerals  to  compressions  of 

V/V  0.7. 
o 

B.  Ultrasonic  Measurements 

Although  some  measurements  of  the  elasticity  of  single  crystal  spinel 
were  carried  out  by  R.  J.  O'Connell  summarized  in  Table  2,  our  major  effort 
was  concentrated  on  a  very  complete  study  of  rock  salt  at  high  temperatures 
and  pressures.  The  work  reported  in  (paper  1])  represents  the  first 
(elasticity)  measurements  of  this  material  at  simultaneously  high  pressure 
and  temperature.  Of  particular  interest  are  the  pressure  dependence  of  the 
derivative  (dK^,/dP)  and  d  K^/dF^  which  were  both  measured  for  the  first  time 
(Fig.  3) .  The  complete  reduction  of  the  ultrasonic  data  in  terms  of  a  p-V-T 
equation  of  state  permitted  calculating  the  variation  of  the  Gruneisen  parameter 
as  a  function  of  volume  and  temperature  Fig,  4.  As  can  be  seen  in  the  figure, 
y  is  nearly  independent  of  temperature,  at  constant  volume,  as  is  assumed, 
in  the  commonly  employed  Mie-Gruneisen  equation  of  state. 
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(dKT/dP) 


Figure  3 


.  (3l\j,/3P)  versus  temperature  at  0,  and  10  kbar  from  the  ultra¬ 

sonic  measurements  on  KaCl  of  Spetzler  et  al.  1972,  paper  11). 
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Figure  4.  Contours  of  the  Gruneisen  function  Y(V,T)  over  the  volume  and 
temperature  range  of  the  ultrasonic  measurements  for  NaCl. 

The  dashed  line  shows  the  zero  pressure  isobar  (after  Spetzler 
ct_al.,  1972,  paper  11). 
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